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Abstract. We introduce and study the block voter model with noise on two- 
dimensional square lattices using Monte Carlo simulations and finite-size scaling 
techniques. The model is defined by an outflow dynamics where a central set of Npcs 
spins, here denoted by persuasive cluster spins (PCS), tries to influence the opinion 
of their neighbouring counterparts. We consider the collective behaviour of the entire 
system with varying PCS size. When Npcs > 2, the system exhibits an order-disorder 
phase transition at a critical noise parameter Qc which is a monotonically increasing 
function of the size of the persuasive cluster. We conclude that how large the PCS is 
more power of persuasion it has. It also seems that the resulting critical behaviour is 
Ising-like independent of the range of the interactions. 
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In nonequilibrium dynamical systems there are usually no energy functions and their 
time evolutions are defined by dynamical rules. These rules can be divided into two 
groups, namely, infiow and outfiow dynamics[lj. For infiow dynamics the center spin 
is infiuenced by its nearest neighbours. A well-known example of such dynamics, in 
nonequilibrium systems, is the majority- vote model in social sciences [21 [3l H]. On the 
contrary, for outfiow dynamics the information fiows from the center spin (or cluster of 
spins) to the neighbourhood. The Sznajd model |5], which was introduced to describe 
opinion formation in social systems, falls on this category. It is based on the fundamental 
social phenomenon called social validation [6]. Figured] illustrates the fiow of infiuence in 
a typical outfiow dynamics where the central spins (circle) try to infiuence the opinion of 
their neighbouring counterparts (rectangles). We shall nominate the central set of spins 
by persuasive cluster spins (PCS), and we aim in this paper to consider the collective 
behaviour of the entire system with varying PCS size. In short, the PCS will try 
to infiuence the neighbourhood and we might expect how large it is more power of 
persuasion it would have. 



n ■ 




Figure 1. Schematic representation of the outflow dynamics showing the central 
cluster of persuasive spins (PCS) and adjacent spins (squares). 

The present block voter model introduces long-range environment-behaviour 
interactions in the system by considering outfiow dynamics. The range of the interaction 
is defined by the number of spins Npcs inside the persuasive cluster. In principle the 
size of PCS can grow until to infinity (infinite range, mean-field limit). However, in the 
limit Npcs — > c>o, in which every spin has the same strength of interaction with each 
other spin on the lattice, the system's behaviour should be predict by the mean-field 
theory with classical critical exponents, namely, /3 = 1/2, 7 = 1, and z/ = 1/2 [7l[8l|9]. 
Here we will not consider this limit and restrict the simulations to finite systems taking 
into account the relative sizes of PCS and the whole system, which limit the number of 
spins in the persuasive cluster. 

The remainder of the paper is organized in the following way: In Sec. 2 we introduce 
the model and describe the methodology used in the simulations. In Sec. 3 we present a 
discussion of our numerical results along with a finite-size scaling analysis of the relevant 



Block Voter Model 



3 



quantities, and we conclude in Sec. 4. 



2. Block Voter Model 



The model system consists of a set of two-state spin variables cxj, associated with 
the zth vertex of a regular square lattice of linear size L = \fN . Each spin can 
have two possible values o"j = ±1 corresponding to the two opposite opinions in a 
referendum. The system's evolution starts from a completely ordered state with all 
spins point to the same direction. At each time step, a block of spins randomly chosen 
tries to influence the spins in the neighbourhood passing to them the opinion of the 
majority of its components. Moreover, all spins have some resistance for accepting such 
outflow influence. So, independently, each spin adjacent of the PCS agrees with the 
PCS majority state with probability Pagree and adopt the opposite state with probability 
q = 1 — Pagree- In terms of the noise parameter q, the rate of flipping for each adjacent 
spin is given by 



wiai 



1-fl 



5=1 



(1) 



where the summation is over all Npcs sites that make up the persuasive cluster, and 
S{x) = sgn{x) if a; 7^ and S{0) = 0. From equation [H we notice that the block voter 
model shares features with the majority-vote model [3111]; e.g., both models exhibit up- 
down symmetry and are endowed with spontaneous broken symmetry on the parameter 
q, having then an order- disorder phase transition. However the majority-vote model is 
a kind of inflow dynamics which takes into account only the opinion of the neighbouring 
spins of a selected node i. Here we consider the majority opinion of a block of spins of 
size Npcs with influence upon its neighbourhood. The block model is a kind of outflow 
dynamics defined in terms of two control parameters: q and Npcs- 

To study the effect of the noise parameter q and the range of interaction Npcs 
on the phase diagram and critical behaviour of the block voter model, we consider the 
magnetization Ml, the susceptibility xl, and the Binder's fourth-order cumulant ?7r p^. 
which are defined by: 



(2) 




sample 



(3) 
(4) 



'time I sample 

where = is the number of spins in the system. The symbols < ■ ■ ■ >time and 
< ■ ■ • >sampie, respectively, denote time averages taken in the stationary state and 
configurational averages taken over several samples. 
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We have performed Monte Carlo simulations on regular square lattices of sizes 
L = 100, 120, 130, 140, 150, 160, and 180, for values of Npcs satisfying the following 
constraint between the size of the persuasive cluster and the whole system: ^^jf^ < 1%. 
Time is measured in Monte Carlo step (MCS), and 1 MCS corresponds to attempts 
of changing the states of the spins. So, for a given PCS we consider two procedures: 
For asynchronous update we choose randomly a spin adjacent to the PCS and try to flip 
it with the probability given by equation ([1]). By repeating this procedure N times, we 
have accomplished one MCS. In the case of synchronous update, however, considering 
that for a given Npcs we have Nadj adjacent spins, we can update simultaneously all 
theses spins with the same probability [equation ([T])]. One MCS is accomplished after 
repeating this procedure = times. For each kind of update, we waited 10^ MCS 
to make the system reach the steady state, and the time averages were estimated from 
the next 10^ MCS. In the critical region larger runs were needed, and we used 3 x 10 
MCS to reach the steady state and 20 x 10^ MCS for time averages. For all set of 
parameters, at least 100 independent runs (samples) were considered in the calculation 
of the conflgurational averages. 

3. Results and Discussion 

Figure [2] shows the dependence of the order parameter and the susceptibility xl 
on the noise parameter, obtained from Monte Carlo simulations on square lattices 
with L = 140 (A^ = 19600) and several sizes of the persuasive cluster spin, namely, 
Npsc = 4,9,16,25,36,49,64. In Fig. 2(a) each curve for Ml, for a given value of L 
and Npcs, suggests that there exists a phase transition from an ordered to a disordered 
state, characterized by a spontaneous broken symmetry at a particular value of the noise 
parameter, namely, q = Qc- In the thermodynamic limit (A^ — t- oo), we should expect 
that below the critical noise Qc the system has a nonzero magnetization, whereas the 
magnetization vanishes for q > qc. For flnite systems, the value of q where each curve 
for Xl in Fig- 2(b) has a maximum is identifled as q^L) for the corresponding Npcs- 
We note that the transition occurs at a value of the critical noise parameter which is an 
increasing function of the size of the persuasive cluster. 

For each Npcs, we can obtain the critical value qc by calculating the Binder's 
fourth-order magnetization cumulant Uiiq) (Equation HI) as a function of the noise 
parameter q, for several lattice sizes L. For sufficiently large systems, these curves 
intercept each other in a single point U*{qc). The value of q where occurs the 
intersection equals the critical noise qc, which is not biased by any assumption about 
critical exponents since, by construction, the Binder's cumulant presents zero anomalous 
dimension |10[ flT]. In Figure [3] we plot the reduced fourth-order Binder's cumulant for 
lattice sizes L = 100, 120, 140, 160, and 180 and four different values of Npcs- From 
each set of curves, we obtain the critical noise parameter qc as well as the critical 
value U*{qc). As we can notice, there exists a strong dependence between the critical 
noise and the size of the persuasive cluster, i.e., when Npcs increases the critical noise 
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Figure 2. Magnetization and susceptibility as functions of the noise parameter q, for 
L = 140 and values of Npcs = 4, 9, 16, 25, 36, 49, 64 (from left to right). 




0,4 



q 



Figure 3. A set of fourth-order reduced Binder's cumulant as a function of q for 
Npcs — 4, 9, 16, 25 (from left to right) and five values of lattice sizes L {L = 
100, 120, 140, 160, and 180). The critical value U* = 0.606 ± 0.004 is shown by the 
horizontal line. 

also increases. Despite the observed dependence between qc and Npcs, the value of 
the Binder's cumulant at the intersection point U* does not depend on the size of the 
persuasive cluster: We obtained U* = 0.606 ± 0.004 (discontinuous horizontal line in 
Figure [3]), for all Npcs, which is in agreement with the result for the Ising model on the 
regular square lattice [I2]. In order to construct the phase diagram for the block voter 
model, we have performed this analysis for several values of the parameter Npcs- 

The phase diagram in the Npcs x Q parameter space is shown in Figur^H The 
curves were obtained using both the synchronous update with all adjacent spins to the 
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Figure 4. The phase diagram of the block voter raodel, showing the dependence of the 
critical noise parameter qc on the size of persuasive cluster Npcsi obtained from MC 
simulations using both synchronous (open symbols) and asynchronous (closed symbols) 
updates. 



persuasive cluster being updated at the same time, and the asynchronous one when 
just a single spin randomly chosen in the neighbourhood of the PCS is updated. For 
both updates, there exists a phase transition only for values of Npcs > 2. Moreover, 
the critical noise parameter Qc is an increasing function of the number of spins in the 
persuasive cluster. In other words, we can conclude that how large the PCS is more 
power of persuasion it has. 

We turn now to the finite-size scaling (FSS) theory [131 El that allows us to 
extrapolate the information available from finite-system simulations to the relevant one 
in the thermodynamic limit. The critical behaviour of the block voter model is given 
by: 

Mi(g) ~ L-^/^MieL^/"), (5) 

XL{q) ~ L^/'x{eL'/n, (6) 
U^iq) ~ UieL'/^, (7) 

where e = (q — qc) is the distance from the critical noise. The exponents /3/z/, 7/z/, and u 
are, respectively, associated to the decay of the order parameter M^^q), the divergence 
of the susceptibility XL^q), and of the correlation length ~ The universal scaling 
functions M{eL^^'^) and ^(eL^/^) depend only on the scaling variable x = eL^^'^. 

The correlation length exponent u can be obtained from the derivative of the 
Binder's fourth-order cumulant with respect to the noise parameter, namely, U'^{qc)- 
Moreover, we can check our estimates for the critical exponents considering the 
hyperscaling relation derived from the Rushbrooke and Josephson scaling laws, namely: 

2/3/z. + 7/z/ = ci, (8) 
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Figure 5. Log-log plot of the size dependence of M^i^q^ for PCS sizes Npcs = 
4,9,16,25,36,49, and 64 from top to bottom. The slopes of the solid lines, obtained 
from the best linear fits to the data points, yield estimates for the exponent P/v in 
good agreement with the exact value 0.125. 
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Figure 6. The same as in Fig. [5l but for the dependence of XL (9c)- The slopes of the 
straight lines give the exponent •^/v close to the exact value 1.75. 
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Table 1. Results for the critical noise Qc, the critical exponents P/v,^ /v,l/v and 
the effective dimensionality D^f / for the Block Voter Model on regular square lattices 
considering different sizes of persuasive cluster spin Npcs- 



Npcs 












4 


0.0846(2) 


0.130(2) 


1.72(4) 


1.05(2) 


1.98(4) 


9 


0.2080(1) 


0.124(4) 


1.76(3) 


0.98(1) 


2.01(4) 


16 


0.2680(1) 


0.114(2) 


1.73(1) 


0.96(2) 


1.95(2) 


25 


0.3130(4) 


0.132(1) 


1.74(1) 


0.95(2) 


2.00(1) 


36 


0.3406(3) 


0.130(3) 


1.70(3) 


0.96(3) 


1.96(4) 


49 


0.3630(1) 


0.113(3) 


1.70(2) 


0.96(2) 


1.96(3) 


64 


0.3790(2) 


0.118(2) 


1.71(4) 


0.96(3) 


1.95(4) 


Obtained using D^ff 


■ = 2/3/zy ^ 









which is vahd for Euchdean dimension d less than the upper critical dimension d^ [T5] . 

Figures E] and E] show the dependence of the magnetization and susceptibility on 
system size L, at g = qc, when different values of Npcs are considered. From the scaling 
relations (equations E]) and (E]), respectively, for a given Npcs the corresponding straight 
lines in these log- log plots have slopes equal to the exponents /3/z/ and 7/z/. This analysis 
yields exponents very close to the two-dimensional Ising exponents, namely, P/v = 0.125 
and ^jv = 1.75. Similarly, from the slopes of the resulting straight lines in the log-log 
plots of the derivative of the Binder's cumulant with respect to the noise parameter, 
at g = gc, we determine the exponent associate with the correlation length in good 
agreement with the exact value v = 1.0. Finally, we have checked whether the calculated 
exponents satisfy the hyperscaling relation (equation |H]) with -De// = + 7/z/. For 
all PCS sizes considered we obtained -De// = 2, within error bars. Our results are 
summarized in Table I. 

According to Grinstein's conjecture [16], all equilibrium models with up-down 
symmetry on regular lattices have Ising-like critical behaviour. This conjecture has 
also been verified for several nonequilibrium model systems p!71 [T8l [iQl 120] . Moreover, 
by considering the universality of the critical exponents for arbitrarily large but finite 
range of the interactions (see, e.g., [211 [22] and references therein), we actually should not 
expect even such a small dependence in the calculated exponents with the number of 
persuasive spins Npcs- In fact, the critical behaviour of the present nonequilibrium 
Block Voter model and the equilibrium (finite-range interaction) Ising model are 
described by the same critical exponents. 

In order to quantify the above statement, we shall consider a more accurate analysis 
that consists of obtaining the universal curves for the mag netization M(x) = ML{q)Ll^''' 
and for the susceptibility x{^) = XLiq)L"'^'^ with x = {q — qc)L^^'^, which follows from 
equations |5] and El respectively, and represent the data collapsing for fixed Npcs and 
different sizes L of the system. It is worth to note that the universal curves are obtained 
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Figure 7. Data collapsing of the order parameter and of the susceptibility for five 
different values oi L = 100, 120, 140, 160, and 180 with Npcs = 9. The universal 
functions are consistent with Ising exponents: P/v — 0.125, ^ /v — 1.75 and 1/v = 1.0. 

once the correct values for the exponents are used. Figur^ shows the data collapse for 
the case of Npcs = 9 and five different system sizes L = 100, 120, 140, 160, and 180, 
using the following exponents: /3/z/ = 0.125, 7/z/ = 1.75 and l/u = 1.0. A similar 
analysis apphed to other values of Npcs supports the conclusion that the present model 
is in the same universality class of the two-dimensional Ising model. 

4. Conclusion 

In this paper, Monte Carlo simulations and finite-size scaling theory was used to study 
a medium-ranged interactions version of voter model with outflow dynamics. The 
model is defined in terms of two parameters, the size of the persuasive cluster of spins 
(Npcs) and the noise parameter q associated with the resistance that every spin has 
for accepting such outflow influence. Considering both synchronous and asynchronous 
updates, the resulting phase diagram in the Npcs x 1 parameter space indicates that 
the region where there exists an ordered phase increases with increasing range of the 
interactions, meaning that how large the persuasive cluster is more power of persuasion 
it has. The calculated critical exponents for different sizes of PCS support the well- 
known criterion of universality class, stating that medium-range interactions models 
with up-down symmetry exhibit Ising-like critical behaviour for all arbitrarily large but 
finite range of the interactions. 
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